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For coupled-channel resonance scattering we derive a model with a closed form solution for the 
T-matrix that satisfies unitarity and analyticity. The two-channel case is handled explicitly for 
an arbitrary number of resonances. The method focuses on the expansion of the transition matrix 
elements, in known analytical functions. The appropriate hadronic form factors and the related 
energy shifts can be determined from the scattering data. The differences between this method and 
the if-matrix and the Breit-Wigner approximation are illustrated in the case of the Sii resonances 
S11 (1535) and S11 (1650). 
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I. INTRODUCTION 

Analyticity is numerically the hardest problem to keep track of in a scattering problem. Both the singularities 
from on-shell states, which yield imaginary parts and the effects of thresholds require careful analysis. This 
is one of the main hurdles for the development of unitary models for coupled-channel scattering. 
In this short paper a set of consistent approximations that yield unitary and analytical results for coupled 
channel problems is described. It is accompanied by the open source program BUBBLEGUM, which yield 
numerical T-matrix results, together with the corresponding if-matrix approximation, Q Breit-Wigner 
approximation, ^ and the perturbative results. The general purpose of both the paper and the computer 
program is threefold. First, to show that approximations that violate unitarity or analyticity will generally 
deviate considerably ||j from exact results in the case of strong interactions and thresholds close to resonance 
values. Second, to show that the exact results can be derived straightforwardly with a reasonable effort, 
in a manner that can improved be systematically. Third, to show that scattering data do not only contain 
pole information of the resonances but also the non-perturbative matrix elements between asymptotic, or 
scattering states, and resonances. Extraction of these matrix elements from the scattering data might lead 
to better constraints on models for hadrons. 

Furthermore, substantial energy shifts occur through coupling of resonances with asymptotic states, or decay 
channels, which can only be analyzed in a fully analytic approach. Since hadrons are composite particles, they 
have form factors that suppress high-energy processes. These form factors and the energy shifts are related 
and should be treated consistently, which neither the .ftf-matrix approach nor the Breit-Wigner approximation 
can do. Depending on the large energy behavior of the form factor, the shift can vary substantially. In the 
case of suppression at high energy, the shift will change sign as scattering energy increases, from negative to 
positive. 

It is common folkore that bare masses have no meaning. It stems from renormalizable field theory, where 
indeed masses can undergo abritrary shifts through finite renormalization. However, in hadron theory 
constituent quark models and other hadronic models yield the bare states with the bare masses and their 
couplings to the decay channels. In this case a comparison between scattering data and these models should 
incorporate the energy shifts through the non-perturbative coupling of bare states and decay channels. 



II. THEORY 



The simplest process contributing to s-channel scattering is the tree diagram, where the two-particle state j, 
in a particular partial wave, forms a resonance state b which decays into a two-particle state i. The transition 
amplitude is given by: 

- 9ibfi(E)yfc(Ej E _ -L _ _ yfcWfjWgh , (1) 

where gn> and gjb are the coupling constants, p(E) is the phase space for each state, f(E) the form factors, 
Mb the mass of the state b, and E the scattering energy. 
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This lowest order diagram is a poor approximation || to the actual scattering process, especially for hadron 
dynamics. First, the coupling constant is large, therefore, during scattering, the initial state j goes, on 
average, through several intermediate states, before ending up in the final state i. Second, for scattering 
energies E close to the resonance mass M, the interaction is strong due to the small energy denominators, even 
if the coupling is small. Hence, the simple tree diagram, which corresponds to the lowest order perturbative 
expansion, fails here. If the problem would be solved consistently, the resonance will acquire a width, which 
corresponds to the shift of the pole into the complex plane. 

Solving the scattering problem consistently means iterating the elementary scattering processes such that it 
reaches a fixed point. There are many more or less equivalent ways to find the self-consistent solution. In 
this paper we will use the Lipmann-Schwinger equation [§, §]: 

- t(°) j-VtI ) 



^2 T cd I7_ TP I j, ^da , (2) 



ca ■ ^-ca E-E d + ie 
d 

which can be written down in many other forms. In practice, this equation determines the transition 
amplitude T, which is the sum over all possible processes, with different numbers of intermediate states. 
For the case of two coupled channels the matrix Lipmann-Schwinger equation will look like: 

in *ia A = y , x (*Wi + Wi W t n t 12 \ 

t 21 t 22 J + ' ^ o iW 2 + W 2 ) ' \ fai h 2 ) ' [6> 
where 

X -\2-.— W a \9* la 92a \92a\ 2 )) ■ (4) 

The reduced transition amplitudes t are defined as: 

T ij (E) = \i)t ij (E)(j\ , (5) 

with the states \i) given by: 

\i) = JdE \fp~[E)fi(E)\(f)i(E)) . (6) 
\<pi(E)) is the asymptotic state i with energy E. 

The solution for this case of two coupled channels can still be given in a simple closed form: 

11 I - XnWi - X 22 W 2 + det[X]WiW 2 ' U 
t = X 22 - det[X]W! 

22 1 - X n Wi - X 22 W 2 + det[X]WiW 2 ' U 

t „ = ^ (g) 

1 - X n Wx - X 22 W 2 + det[X}WiW 2 ' yj 

21 l-XiiWi-X 22 W2 + det[X]WiW2 ' 1 ' 

where Wj is defined as: Wj = iWj(E) + Wj(E), and the functions W are defined through the following 
integrals: 

Wi(E) = 7T / dE'{<j> i [E')\<t>i{E))p i {E)f i {Ef = -Kp i {E)f i {E) 2 , (11) 



,r,,D = i'v / dE' Pi{E E ] ^' )2 , (12) 



where PV stands for a principal value integral, and the asymptotic states are orthonormal: {<j}i{E')\4>i{E)) = 
5(E' — E). It is clear that these results can easily be extended to incorporate an arbitrary number of channels. 
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Both p and / are real; any phase information is to be incorporated in the coupling constants. In the absence 
of a form factor / the second integral generally does not converge. It requires regularization, and possibly 
a motivation for this regularization through renormalization. Subtracted dispersion integrals effectively 
correspond to regularization through local subtractions. Another option is to set W to zero, and only use 
the phase space for the imaginary part, which is known as the isT-matrix method. 

For a general case the dispersion integral Eq. ( |l2| ) cannot be solved analytically. However, since the form 
factor is rather arbitrary, it can be chosen such that pi(E)fi(E) 2 has a known solution for the dispersion 
integral. An expansion in a set of these functions, allows one to determine p(E)f(E) 2 as part of the analysis 
of scattering data of resonances. This quantity can, and should, be compared to models of hadrons since the 
function p(E)\gf(E)\ 2 is the square of the transition matrix element of the Hamiltonian: 



where M is the bare resonance state. The transition matrix element is expected to fall-off at higher energies, 
due to the composite nature of hadrons. 

Another constraint on W(E) is given by general covariance, which requires that W(E) is analytic in s = E 2 
with all the cuts on the real axis associated with decay channels, or asymptotic states. In Feynman pertur- 
bation theory covariance is the result of summing the different time-ordered diagrams. In the language of 
states, these different time orderings correspond to different states, which often correspond to decay channels 
that will be of little significance in the kinematical region of interest. However, as restoring covariance comes 
at little cost, we replace the energy dispersion relation with: 



would yield generally a small difference in W. This, however, includes the second time ordering and restores 
analyticity in s, which is a consequence of full covariance. Moreover, analyticity in s implies a restriction 
on possible approximations for W(E), i.e., there should not be a cut at s = 0. In W(E) itself the second 
time ordering is not included, as it would correspond to a different decay state which is generally of little 
consequence, however, would yield additional singularities to keep track of. This feature of manifest rela- 
tivistic invariant formulation makes the evaluation of higher order Feynman diagrams in Minkowski space 
complicated. The advantage of the one-state-one-singularity approach advocated here is that the singular 
structure of the perturbative kernel of the Lipmann-Schwinger equation is near to trivial. Characterizing 
states by their energy simplifies the equations to such an extent that a greater part of the calculations can be 
performed analytically, as we will see below. Note, when the kinematical domain of validity is restricted to 
the energies where only the designated states, that appear as channels in the Lipmann-Schwinger equation, 
can go on shell, the result is fully covariant with the replacement Eq. (|l4|). 



Given a particular threshold behavior n, associated with a partial wave, we use as leading order approximation 
for W{E): 



where Sth is the threshold energy squared, and [i] is the integer part of i. The power n is not only related 
to the partial wave, but also to kinematical factors which are different for a heavy-light system like irN and 
an equal mass system, like ttit. Furthermore, it depends as well on the kinematical range of the threshold, 
i.e., if the threshold extends to where the masses are comparable to the energy, it would yield a different 
threshold behavior. For example, in the ttN scattering in the A region, the real threshold behavior stops at 
about 50 MeV above the threshold. 

An additional expansion to model or fit the function W away from threshold is given by a polynomial P m 
of order m in £ = (s — s t h)/s: 



p(E)\gf(E)\ 2 = ME)\H\M}\ 2 



(13) 




(14) 
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(15) 
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i n/2 Pm(0 



(16) 
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(17) 
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0(s/sth ~ l)\/ith y2 c i w {n/2+i)([n/2]+i+l) , (18) 



i=0 



such that W is analytical in s — and has the dimension of energy. The first coefficient Co = 1 such that the 
threshold behavior is determined solely by the power n and the coupling constant. The coupling constants 
gib are dimensionless. Note, that no additional dimcnsionful quantities are introduced beside the threshold 
energy, which makes a good candidate for a normalized description of resonance scattering, which is up to 
now littered with different form factors with all kind of dimensionful quantities with little meaning. If any 
dimensional quantities, such as the QCD scale, pion decay constant, or the pion mass, are relevant the can 
be expressed as a number times the threshold energy ^/sth- 

Each of these terms in the polynomial expansion has a closed form solution for the dispersion integral, where 
we distinquish integer and half- integer values for n in w nk : 

( s — IV 

w nk = 6(s-l)ir {S „ fc > (19) 



*W = (l-To) **- 1 ^' 1 -*' (20) 
Re(s-l)( 2n+1 )/ 2 

w^±l k = f k (21) 

w^ k = (-!)»(! -r )7r^ % (22) 

z s 

where Tq refers to a Taylor-Laurent expansion of all the singular terms in s , and s = s/sth- The function 
£ m /s Eq. ( |l6| ) peaks at s = (m + I)sth, while £, m /V$ peaks at s = (2m + l)s t h, such that every higher order 
term in P m (0 extend the range of energies over which W(E) can be fitted. The polynomials form a L 1 basis 
for W on < 1/(1 + < 1- 

The real part w can be derived in closed form by expanding the series and resumming the terms, analogous 
to the calculation in jlC| |. However, since the singularity at s — has to cancel, and the difference can 
only be meromorphic functions that fall off at infinity, it is simple to see that the real part is restricted to 
the trivial analytical continuation and the Laurent series. The functions w generally have a large negative 
value, leading to a negative mass shift. This is mainly due to the large high-energy tail W 2 n +i k ~ 1/%/* or 
w n k ~ 1/s. This leading order behavior is necessary to be able to fit arbitrary W(E) dependence, however, 
might be suppressed in the actual scattering, with an equal reduction of the mass shift. 
Using these functions to solve the Lipmann-Schwinger equation in closed form effectively means making a 
simultaneous, or coherent, bubble sum approximation in each of the decay channels for each of the reso- 
nances. Therefore we refer to this method as the Coherent Bubble Sum Approximation. Single bubble sum 
approximations have been used widely. It corresponds to dressing the propagator of the resonance with its 
single or multiple decay channels. However, the case of an unitary, analytical, and covariant approach with 
channel mixing has only appeared at the heart of a few coupled-channel analyses with fixed form factors. 
|?], ||, |j| In the case of a single resonance, the bubble sum can be seen as a dressing of the resonance par- 
ticle; a change in the particle properties. In the case of coupled channels such an interpretation does not 
exist. In this case only the physical observables in the scattering experiment are free of ambiguities in the 
interpretation. In simple terms, the resonance properties are intertwined and co-dependent. 
Solving the Lipmann-Schwinger equation Eq. (|^) reduces to solving the algebraic equation Eq. (|3[), and 
determining the Laurent-Taylor expansions, which is automated in the BUBBLEGUM code. 



III. MASS SHIFTS 



The if-matrix method leaves the real part of the resonance pole in the same position. The coherent bubble 
sum shifts the pole in order to satisfy analyticity contraints. These shifts can be substantial. However, the 
word "mass shifts" suggest that all effects can be reabsorbed in shifting the mass such that it accounts for 
main differences between if-matrix and T-matrix results. This is generally not the case. In a coupled channel 
problem, the mass shifts differ from channel to channel. In the if-matrix the real part of the pole, i.e., the 
position where the real part of the amplitude crosses the axis, is at the same location in every channel as 
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long as the coupling constant is not too large. Only for large coupling constants unitary yields stringent 
constraints for the ii"-matrix, which can shift theses locations independently. For the exact result this does 
not hold. The isobar model, jj| which implies resonance properties independent of the decay channels, is 
therefore flawed from the start. The procedure of finite mass renormalization is not unique. 
However, in order to make the comparison between if-matrix and exact results more meaningful for the 
if-matrix method, one can make a global mass shift AM a for resonance a, which, in some way, averages the 
shifts in each of the channel: 

A ,, r {E-M a )(l-X u W 1 -X 22 W 2 + det[X]W 1 W 2 ) 

AAi„ = lim — * , (26) 

E^M a ^§_(E - M a ){\ - XuWi - X 22 W 2 + det[X]WiW 2 ) 

where the expression is given by the common denominator in all of the T-matrix channels, Eq. (|7-10). 
The term in the denominator corresponds to the appropriate wave function renormalization. In the case of 
small coupled-channel effects and small overlap of the resonances the real part of the shifted pole positions 
coincide with the expected channel-independent pole position. However, in the case large couplings or 
coupled-channel effects this independent resonance mass shift breaks down, i.e., linear approximation of the 
denominator around the pole is no longer valid. In that case an iterative optimizing scheme to recover the 
mass shifts is required. However, as the mass shift is not an invariant quantity, the results from such a 
scheme would have no real significance. Eventually one has to accept that the scattering data and the bare 
properties and couplings are the only well-defined quantities. However, the bare masses must always be 
given in combination with the coupling constants and the form factors. 



IV. t-EXHANGE WITHOUT OPEN CHANNELS 



This method of comparing scattering data to resonance properties can also serve as an approximate method 
for handling additional i-exchange processes. For low and medium energy scattering, i-exchange diagrams 
will not go on-shell, i.e., there is no imaginary part associated with the perturbative diagram T c ^ h . However, 
the real part is a direct consequence of the energies at which the t-exchange diagram goes on-shell and strength 
in the channel. Therefore, the contributions can be approximated by a single resonance: 

ImT ( x °i = - J ^E')dE'g ix f i (E,E')y/^nS(E-E')y/^f j (E,E')g^ 

-mxfi(E, M cB )y/ Pi (E)n6(E - M cS )^ Pj (E)f 3 (E, M eS )g* x , (24) 

T o°L = miJi (E, M eS ) Jp~(E) 1 Jp~{E)f 3 (E, M cS )g* x , (25) 

£j — -M c ff + v 

where [J-(E') is a measure of the kinematical range over which the i-exchange diagram can go on shell, and 
M e ff is the mean of that measure: 

J (i(E')dE' = fx , (26) 

J n(E')dE'E' = M cff . (27) 

In the case where the t-exchange diagram has an imaginary part in the kinematical range of interest, this 
approximation will not hold. More appropriate methods for this problem are under investigation and will 
be reported later. 

In effective field theory highly virtual intermediate states are replaced by effective contact interactions. The 
effective resonance described here does more than that. The effective resonance that replaces the virtual 
i-exchange incorporates also the leading energy dependence of the i-exchange contribution. 



V. Sn RESONANCES 



The 5ii resonances at 1.535 GeV and 1.65 GeV are a good test case for the coupled channel analysis of 
overlapping resonances as the threshold of the rjN channel is in the same region. However, there is only 



G 




1.25 1.5 1.75 1.25 1.5 1.75 

Energy [GeV] Energy [GeV] 



FIG. 1: The transition amplitudes in Sn nN scattering. The top row are the real and imaginary part of the nN — » 7riV 
amplitudes, and the bottom row the ttN — > rjN amplitudes. The dots and circles are partial wave data, without error 
bars, the solid lines are the T-matrix results, the dashed lines the if-matrix approximation of these results, with 
identical parameters, and the dot-dashed lines the Breit-Wigner results, as described in the paper. 



limited data available for the rjN channel. We fit the data between threshold and 1.77 GeV with minimal 
model of two bare masses and four coupling constants. Beyond the resonance region there is still a substantial 
amplitude, and little data in other than the irN channel, which might indicate the importance of ttttN final 
states, and f-channel exchanges. A kink in the ttN — » ttN data at 1.3 GeV is a further indication of ttttN 
states. For the clarity of the example we do not include the high-energy region and work with the minimal 
parameter set of six parameters. We use the GWU-VPI partial wave datajLlJ for the ttN elastic amplitude 
and the results from the Pitt-ANL analysis || for the ttN —>■ rjN channel. Standard fitting algorithms do 
not work properly as the variables are highly correlated in a non-trivial manner due to the unitarity and 
analyticity conditions. Instead we performed a global search, with an increasing mesh, zooming in on the 
optimal fit. 

In Figure [l] we compare the exact T-matrix results with the if -matrix and Breit-Wigner approximations. 
The top left figure shows real part of the irN — > ttN transition amplitude, while the imaginary part appears 
on the right. The bottom row is the same combination for the irN — > rjN transition amplitude. The bare 
masses are 1.60 GeV and 1.79 GeV, while the renormalized masses come to 1.51 GeV and 1.71 GeV. The 
latter are used for the if -matrix and the Breit-Wigner approximations. Without these mass shifts Eq. (p3|), 
the approximate results would have no resemblance with the data. The dotted line and the circles are the 
Pitt-ANL and GWU-VPI partial waves for the rjN and the irN channels respectively. The deviation for 
the rjN partial wave is due to the large error uncertainty, which led to the dominance of the 7riV channel 
in the fit. However, varying the relative importance of the rjN data, with respect to the 7r/V data, did not 
significantly altered the results. It seems that the Pitt-ANL and GWU-VPI partial wave analysis, which 
include background terms, lead to an inconsistent rjN amplitude that cannot be fitted with with a resonance- 
only model like the coherent bubble sum approximation. With the current rjN data it is not possible to 
determine the nature of the deviations, and whether the data favors a background contribution. However, 
there seems to be no natural explanation for a background contribution. The background should model some 
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degrees of freedom, or states, that are important yet virtual, i.e., not associated with open channels, such 
as, perhaps, the p-exchange between the irN pair. 

Including the first term c\ in the expansion Eq. ( |l8| ) of the form factors halves the mean square deviation 
between the data and the fit. Notably, it increases the effective width of the irN form factor, while it 
decreases the effective width of the rjN form factors, corresponding, respectively, to a positive and a negative 
expansion coefficient c\. 

Note that it would be possible to obtain better fits in the iCmatrix approximation and the Breit-Wigner 
approximation. However, the purpose of this paper is to point out the differences which are solely the result 
of the approximations made, keeping the model identical. In this case it is clear that the first resonance 
Sii(1535) is overestimated in both approximations. The coupling constants are gi^N = —0.11, g2,-nN — 
0.175, gi, vN = -0.10, and g 2 , V N = -0.11. 



The if-matrix approximation does not require the solution of the dispersion integral. Therefore, it is often 
used to incorporate more complex interactions. In the Coherent Bubble Sum Approximation it correspond 
to setting Wi to zero. The Breit-Wigner approximation is generally not a well-defined procedure. Eventually 
it boils down to replacing the complex interaction by a sum of Breit-Wigner forms with partial widths for 
each of the decay channels in each of the resonances: 



where, in principle, each of the coupling constants and p^s are unknown and should be recovered from the 
analysis. However, given the W(E)'s and the coupling constants, the closest general analogy is given by 
T(E) = Wi(E). From the example of the Sn resonances, where the coupling strength is moderate, it is clear 
that both the if-matrix and the Breit-Wigner approximation leads to significant deviations in the transition 
amplitudes. 

A feature of scattering data where approximate methods such as Breit-Wigner and X-matrix approximations 
yield erroneous results is the interaction between thresholds and subthreshold resonances. A resonance with 
a strong coupling to a channel with a threshold at a higher energy might lead to anomalous threshold 
behavior. Once a resonance has an energy close enough to the threshold energy, it might get attracted into 
the channel, giving either a full circle in the Argand plot, or a sharp spike in the real part of the transition 
amplitude. Neither feature is properly reproduced in the Breit-Wigner or the if-matrix method. 
In this paper the form of W(E) is not derived. Only simple assumptions about its threshold behavior are 
made. The energy dependence can be modeled by a polynomial in the variable (s — s t h)/s. This freedom 
is an essential part of the Coherent Bubble Sum Approximation. If possible, the data should determine 
the function W(E). Its value should be extracted from the scattering data without model dependence. 
It was only assumed that W(E) falls off at infinity, which is considered realistic for composite hadronic 
systems. Analyzing resonance scattering data and modeling transition matrix elements of hadronic states 
are two separate problems and should be treated as such. In future studies we will examine ^-exchange, 
multi-particle final states, and field-theoretical approaches which require renormalization. 
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